Second and higher harmonics generation with memristive systems 
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We show that memristive systems can be used very efficiently to generate passively both double 
and higher frequency harmonics. A technique for maximizing the power conversion efficiency into 
any given harmonic is developed and applied to a single memristive system and memristive bridge 
circuits. We find much higher rates of power conversion compared to the standard diode bridge, with 
the memristive bridge more efficient for second and higher harmonics generation compared to the 
single memristive system. The memristive bridge circuit optimized for second harmonic generation 
behaves as a two-quarter-wave rectifier. 



In nonlinear optics, the phenomenon of second- 
harmonic generation (SHG), demonstrated in 1961 
refers to the possibility of creating an outgoing light beam 
with double the frequency of the impinging bearn^. This 
is usually achieved with the help of nonlinear crystals. 
The smallness of the nonlinearity, however, limits the 
efficiency of direct SHG. The use of lenses and mirrors 
improves the efficiency up to 85% by making the light 
pass repeatedly through a nonlinear medium 3 . A com- 
mon application of SHG in optics is the creation of laser 
sources with modified frequencies. 

In the field of electronics, SHG is termed as frequency 
doubling. Electronic frequency doublers can be divided 
into active ones, requiring an external power source, and 
passive ones, operating exclusively on the power of the 
input signal (see, e.g., Ref. 0). A straightforward re- 
alization of a passive frequency doubler is based on a 
diode bridge^, which produces a full-wave rectified out- 
put (mathematically, the absolute value of a sinusoidal 
input). Using Fourier analysis, it can be shown that the 
diode bridge transforms 4.50% of the input power into 
the second harmonic, and 18.9% of the input power into 
all higher harmonics^. In addition, by symmetry its out- 
put signal contains only even harmonics. 

Memristors 6 , which belong to the more general class of 
memristive systems^, are resistive circuit elements with 
memory. A recent understanding 8 of the memristive 
nature of resistance switching in memory cells has at- 
tracted a lot of attention to the field of memory ele- 
ments 9 , which comprises also memcapacitive and rnemin- 
ductive systems^, namely, capacitors and inductors with 
memory. 

In this paper we demonstrate that a passive circuit 
with a single memristive device generates second and 
higher harmonics signals on a load resistor with signif- 
icantly higher efficiency than that of the diode bridge. 
Moreover, the efficiency of harmonics generation can be 
further improved employing a memristive bridge intro- 
duced in this paper. We note that SHG in large memris- 
tor networks was anticipated in Ref. |Tl| but no detailed 
analysis was given. 

We start by showing that memristive systems indeed 
generate second and higher harmonics using a simple 



model of resistance switching introduced in Ref. |8| to 
describe the response of a thin Ti02 film of a width D 
sandwiched between two Pt electrodes. The phenomenon 
of resistance switching in this system is understood^ as 
a field-induced migration of oxygen vacancies changing 
the width w of a doped low-resistance region. Electri- 
cal current in the positive direction (towards the thick 
line in the circuit symbol of the memristor in Fig. QJa)) 
increases w, decreasing the resistance. If the maximum 
and minimum resistance values, obtained at w = and 
w = D, are denoted by R g and R on , respectively, and 
the charge flow through the device needed to completely 
switch it from one limiting state into another is qo, then 
the device memristance (memory resistance) can be writ- 
ten as& 

w 

R{w) = R oB + (R on - R oB ) -, (1) 

where w/D = q/qo, and q is the cumulative charge flow- 
ing through the device. 

We first consider a sine voltage source of period T con- 
nected in series to a memristive device R(w) and a load 
resistor Ri (Fig. [U^a)). Kirchhoff's voltage law for this 
circuit is 

V sin LJt-q[R(w)+R 1 ] = 0, (2) 
where R(w) is given by Eq. (JTJ and to = 2ir/T. Solving 




FIG. 1. (a) Circuit consisting of a time-dependent voltage 
source connected to a memristive system with memristance 
R(w) and a load resistor R\. (b) Memristive bridge circuit. 
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where I(t) = g, and i?2 = Vo/(gow) is a constant with di- 
mensions of resistance. First, we see that the current I(t) 
in Eq. has the same periodicity as the source and, 
therefore, can be written as a Fourier series. Second, the 
numerator in Eq. (J3J) is a static resistivity current ex- 
pression, while the denominator contains corrections due 
to memory. These corrections give rise to all higher har- 
monics. In many memristive devices R on <C i? ff, and 
it is reasonable to assume R\ < R g. Under such as- 
sumptions, it follows from Eq. ([3]) that the condition 
for considerable generation of higher harmonics is hav- 
ing i? 2 / i? ff of order one (while keeping the denominator 
real). For a given memristor, the ratio R2/R0S can be in- 
creased either by increasing the amplitude of the source, 
Vo, or by decreasing its frequency, lu/2it. 

The memristor model given by Eq. ([1]) combined with 
w/D = q/qo is convenient for analytical calculations. 
However, it does not limit R(w) between R on and R s- 
In order to obtain quantitative results, we then suggest 
a more realistic model consisting of Eq. ([1]) and 

w = — [6(1 - w/D)9{q) + 6{w/D)0{-q)] q, (4) 

where the 9 functions in Eq. (UJ constrain w to satisfy 
< w < D, in agreement with experimental data 13 ' 14 . 

It follows from Eqs. (01 and ([T]) that w has the pe- 
riodicity of q, and that the current I{t) has the same 
periodicity as that of the voltage source. However, I(t) 
cannot be found in a closed analytical form now, and 
we solve the problem numerically. If we denote the av- 
erage power dissipated on the load resistor in the i-th 
harmonic as Pi and the average power of the source as 
^source, three important ratios come to mind: P2 / 'Psource) 
-P3/-P S ource and (J2kLi P k)/ Psource, which are the second, 
third and higher harmonics conversion efficiencies, re- 
spectively. Identifying the system parameters maximiz- 
ing these ratios will enable us to generate, upon appli- 
cation of a suitable band pass filter at the output, the 
desired harmonics passively and with minimal losses. 

Considering all possible parameters, we can write the 
general functional dependence of the above ratios as 



P k 



^source 



= f k (u,Vo,q ,w(t = Q),R ot f,R on ,Ri) , (5) 
= g(u,V ,qo,w(t = 0), R oS , R on , i?i) . (6) 



As noted above, these ratios cannot be written in closed 
analytical form. Moreover, the numerical maximization 
is difficult due to the large number of parameters. It is 
shown below that it is reasonable to select w(t = 0) = 0. 
The number of parameters can be further reduced utiliz- 
ing the Buckingham ir theorem of dimensional analysis 1 ^. 
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FIG. 2. (color online) Current I(t) and corresponding power 
ratio distribution (Pi /-Psource) after maximization of the power 
conversion into the second (a,b), third (c,d), and all higher 
harmonics (e,f) in the single memristive device circuit and into 
the second harmonic in the memristive bridge circuit (g,h). 
I(t) has been obtained using Vb = IV, qo = 3 ■ 10 _12 C and 
R oS = 2(M1 

Its application results in functions of dimensionless vari- 
ables for the above ratios: 
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The left-hand sides of Eqs. ([7]) and © are maximized 
numerically considering the following realistic ranges of 
dimensionless parameters®: 10~ 5 < i?i / i? ff < 10, 
1CT 4 < R2/R0S < 10 and 10 < R oS /R n < 10 3 . We 
find that the maxima of P2/ -Psource and P3/P SO urcc de- 
pend weakly on R s/R n and lie in the range 150 < 
Ros/Ron < 1000. Consequently, we set R s/R on = 200. 
For J2T=2 p k /-Psource it is found that R of f/R on = 1000 
gives the best value. 

Table Q] presents conversion efficiencies and lists the 
optimal parameter values. Clearly, in all cases the sin- 
gle memristive device circuit provides significantly higher 
conversion rates than the diode bridge. Figs. Hla)-(f) 
show the current through the load resistor and corre- 
sponding harmonics power spectrum. Since the memris- 
tive device is a dynamically adaptive system, I(t) is very 
different from the sine shape of the source. 
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TABLE I. Optimized parameter values for single memristive 
device circuit, power conversion rates for diode bridge (Rate, 
4 diodes), optimized power conversion rates for single mem- 
ristive device (Rate, 1 memr.) and memristive bridge (Rate, 
4 memr.) circuits. 



It is also important to check how the initial value 
of w(t = 0) affects the higher harmonics generation. 
We have performed extensive numerical simulations and 
found that the power conversion rates (Eqs. (J5J) and ©) 
always increase with a decrease of w. This observation 
justifies our choice of w(t — 0) = corresponding to 
R(w(t = 0)) = R ff- Moreover, we note that once a 
certain R(w(t — 0)) is selected, the value of Roff be- 
comes not important if R(w(t)) > R on . In order to prove 
this statement, we notice that the load current has the 
same period as the source, which implies that the mem- 
ristance also has such periodicity. The memristance in- 
creases/decreases with negative/positive current, which 
changes sign only at t = nT and (n + 1/2)T, where n is 
an integer. Since the memristance decreases immediately 
after t = 0, we conclude the memristance has maxima at 
t = nT , minima at (n + 1/2)T and no other extrema. 
Therefore, R on < R(w) < R(w — 0) during the entire 
period, implying that having wi = w(t = 0) > and 
off resistance Roff is equivalent to having w(t = 0) = 
and off resistance R(wi). The decrease in the memris- 
tance range reduces the circuit nonlinearity, translating 
into lower efficiency of higher harmonics generation. 

The efficiency of the second and higher harmonics gen- 
eration can be increased even further by utilizing a mem- 
ristive bridge having the geometry of the diode bridge 
(Fig. [lib)). The memristive bridge consists of four iden- 
tical memristive devices described by Eqs. ([I) and (|4| 
rectifying the input signal via the delayed switching ef- 
fect. This rectification reduces the load current peri- 
odicity to T/2, leaving only the even harmonics. Conse- 



quently, the second and higher harmonic generation effi- 
ciencies are improved relative to those of the single device 
circuit, while third harmonic generation is excluded. 

In complete analogy with the single device circuit, di- 
mensional analysis for the memristive bridge leads to 
Eqs. ([7]) and © for the efficiency of harmonics gener- 
ation. Dimensionless parameters maximizing f^/f source 
are R 2 /R « = 2.82, Ri/R oS = 0.0267 and R f f /R on = 
1000 giving P2/P S omco = 40.3%, a substantial improve- 
ment over the single device circuit, with same initial con- 
ditions. I(t) and the power distribution of harmonics for 
the memristive bridge are given in Figs. [2Ig,h). It is 
found that with an increase in R ff/R on , the optimized 
current shape tends towards | sinwij in even quarters of 
each period and zero in odd ones. Such a signal has 
f2/-Psourcc = 45.0%, which is indeed the asymptotic effi- 
ciency at Roff /Ron — > 00 ■ We recognize this signal shape 
as being two- quarter-wave rectified. For higher harmon- 
ics generation, we find the value after optimization of 
(J2T=2 -Pfc ) /^source to be 56.4%, a modest improvement 
over the single device circuit. 

In conclusion, we have demonstrated the potential of 
memristive systems for passive second and higher har- 
monics generation. Using an approach to maximize the 
rate of power conversion for a specific harmonic, we have 
shown that memristive circuits are much more efficient 
for harmonic generation purposes (at optimal operation 
conditions) than the traditional diode bridge. In addi- 
tion, the operation voltage in memristive circuits can 
be lower than that used in diode circuits because of 
0.7V barrier voltage of silicon p-n junctions 5 . We also 
anticipate that memristive devices can be beneficially 
used in harmonics generation in combination with ac- 
tive circuits. An example of such a circuit is a memris- 
tive bridge operating with a high resistance load followed 
by an operational amplifier. Finally, memcapacitive and 
meminductive systems^^ can be used for passive (and 
low-dissipative) higher harmonics generation instead of 
memristive ones. The results of our investigation can be 
readily tested experimentally. 
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